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We propose a classical solution for the kinetic description of matter falling into a black hole, 
which permits to evaluate both the kinetic entropy and the entropy production rate of classical 
infalling matter at the event horizon. The formulation is based on a relativistic kinetic description 
for classical particles in the presence of an event horizon. An H-theorem is established which holds 
for arbitrary models of black holes and is valid also in the presence of contracting event horizons. 
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The remarkable mathematical analogy between the laws 
of thermodynamics and black hole (BH) physics following 
from classical general relativity still escapes a complete 
and satisfactory interpretation. In particular it is not 
yet clear whether this analogy is merely formal or leads 
to an actual identification of physical quantities belong- 
ing to apparently unrelated framework. The analogous 
quantities are E <-> M, T aK and S <-> {l/8na)A, 
where A and k are the area and the surface gravity of 
the BH, while a is a constant. A immediate hint to be- 
lieve in the thermodynamical nature of BH comes from 
the first analogy which actually regards a unique physical 
quantity: the total energy. However, at the classical level 
there are obstacles to interpret the surface gravity as the 
BH temperature since a perfectly absorbing medium, dis- 
crete or continuum, which is by definition unable to emit 
anything, cannot have a temperature different from ab- 
solute zero. An reconciliation was partially achieved by 
m 1975 by Hawking I"], who showed, in terms of quan- 
tum particle pairs nucleation, the existence of a thermal 
flux of radiation emitted from the BH with a black body 
spectrum at temperature T — hK/2TTkB {Hawking BH ra- 
diation model) . The last analogy results the most intrigu- 
ing, since the area A should essentially be the logarithm 
of the number of microscopic states compatible with the 
observed macroscopic state of the BH, if we identify it 
with the Boltzmann definition. In such a context, a fur- 
ther complication arise when one strictly refers to the 
internal microstates of the BH, since for the infinite red 
shift they are inaccessible to an external observer. An ad- 
ditional difficulty with the identification S (l/87ra)^, 
however, follows from the BH radiation model, since it 
predicts the existence of contracting BH for which the 
radius of the BH may actually decrease. To resolve this 
difficulty a modified constitutive equation for the entropy 
was postulated |^,|^, in order to include the contribution 
of the matter in the BH exterior, by setting 



(5' denoting the so-called Bekenstein entropy) where S 
is the entropy carried by the matter outside the BH and 

Sbh 

consequence a generalized second law 



3^T?F identifies the contribution of the BH. As a 



6S' > 
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S' = S- 
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was proposed 0, 1^ which can be viewed as nothing more 
than the ordinary second law of thermodynamics applied 
to a system containing a BH. From this point of view one 
notices that, by assumption and in contrast to the first 
term S, Sbh cannot be interpreted, in a proper sense, as a 
physical entropy of the BH, since, as indicated above, it 
may decrease in time. However, the precise definition and 
underlying statistical basis both for S and Sbh remain 
obscure. Thus a fundamental problem still appears their 
precise estimates based on suitable microscopic models. 
Since the evaluation of Sbh requires the knowledge of 
the internal structure of the event horizon (excluding for 
causality the BH interior), the issue can be resolved only 
in the context of a consistent formulation of quantum 
theory of gravitation This can be based, for exam- 

ple, on string theory [3 and can be conveniently tested 
in the framework of semiclassical gravity J, .8| • Regard- 
ing, instead the entropy produced by external matter S, 
its evaluation depends on the nature, not yet known, of 
the BH. However, even if one regards the BH as a purely 
classical object surrounded by a suitably large number 
of classical particles its estimate should be achievable in 
the context of classical statistical mechanics. 

In statistical mechanics the "disorder" characterizing a 
physical system, classical or quantal, endowed by a large 
number of permissible microstates, is sometimes conven- 
tionally measured in terms of the so-called Boltzmann 
entropy Sb = K\nW. Here K is the Boltzmann con- 
stant while is a suitable real number to be identified 
with the total number of microscopic complexions com- 
patible with the macroscopic state of the system, a num- 
ber which generally depends on the specific micromodel 
of the system. Therefore, paradoxically, the concept of 
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Boltzmann entropy does not rely on a true statistical 
description of physical systems, but only on the classifi- 
cation of the internal microstates (quantal or classical). 
As is well knownQ, Sb can be axiomatically defined, de- 
manding (i) that it results a monotonic increasing func- 
tion of W and (ii) that it satisfies the entropy additivity 
law Sb{WiW2) = Sb{Wi)+Sb{W2). Boltzmann entropy 
plays a crucial role in thermodynamics where (i) and (ii) 
have their corresponding laws in the entropy nondecreas- 
ing monotonicity and additivity. Since in statistical me- 
chanics of finite system it is impossible to satisfy both 
laws exactly, the definition of Sb is actually conditioned 
by the requirement of considering systems with W » 1 
(large physical systems). 

An alternate definition of entropy in statistical me- 
chanics is the one given by the Gibbs entropy, in turn 
related to the concept of Shannon information entropy. 
In contrast to the Boltzmann entropy, this is based on a 
statistical description of physical systems and is defined 
in terms of the probability distribution of the observable 
microstates of the system. In many cases it is sufficient 
for this purpose to formulate a kinetic description, and 
a corresponding kinetic entropy, both based on the one- 
particle kinetic distribution function. In particular, this 
is the case of classical many-particle systems, consist- 
ing of weakly interacting ultra relativistic point particles, 
such as those which may characterize the distribution of 
matter in the immediate vicinity of the BH exterior. 

The goal of the paper is to to provide an explicit 
expression for the contribution S, which characterizes 
Bekenstein law (^, to be evaluated in terms of a suit- 
able kinetic entropy, and to estimate the corresponding 
entropy production rate due to infalling matter at the 
BH event horizon. In addition we intend to establish an 
H-theorem for the kinetic entropy which holds, in prin- 
ciple, for a classical BH characterized by event horizons 
of arbitrary shape and size and even in the presence of 
BH implosions or slow contractions. This is obtained in 
the framework where the classical description of outside 
matter and space is a good approximation to the under- 
lying physics. The basic assumption is that the matter 
falling into the BH is formed by a system of TV ^ 1 clas- 
sical point particles moving in a classical spacetime. We 
adopt a covariant kinetic formalism taking into account 
the presence of an event horizon and assuming Hamil- 
tonian dynamics for the point particles. The evolution 
of such a system is well known and results uniquely de- 
termined by the classical equations of motion, defined 
with respect to an arbitrary observer O. To this purpose 
let us choose O, without loss of generality, in a region 
where space time is (asymptotically) flat, endowing with 
the proper time r, with r assumed to span the set ICR 
(observer's time axis). Each point particle is described by 
the canonical state x spanning the 8— dimensional phase 
space r, where x = {r^,p^). Moreover, its evolution is 
prescribed in terms of a suitable relativistic Hamiltonian 



H = -ff (x), so that the canonical state x = {r^^,p^) re- 
sults parameterized in terms of the world line arc length 
s (see 10]). As a consequence, requiring that s — s{t) 
results a strictly monotonic function it follows that, the 
particle state can be also parameterized in terms of the 
observer's time r. To obtain the a kinetic description for 
a relativistic classical system of N point particles we in- 
troduce the kinetic distribution function for the observer 
O, /3g(x), defined as follows 

Pg(x) = p{x)S{s - s{t))S{^^ - 1) (3) 

where p (x) is the conventional kinetic distribution func- 
tion in the 8— dimensional phase space. Notice that the 
Dirac deltas here introduced must be intended as phys- 
ical realizahility equations. In particular the condition 
placed on the arc length s implies that the particle of 
the system is parameterized with respect to s(t), i.e., it 
results functionally dependent on the proper time of the 
observer; instead the constraints placed on 4- velocity im- 
plies that is a tangent vector to a timelike geodesic. 
The event horizon of a classical BH is defined by the 
hypersurface ru specified by the equation 

R{x) = ru (4) 

where x denotes a point of the space time manifold, while 
R{x) reduces to the radial coordinate in the spherically 
symmetric case. According to a classical point of view, 
let us now assume that the particles are "captured" by 
the BH (i.e., for example, they effectively disappear for 
the observer since their signals are red shifted in such a 
way that they cannot be anymore detected 113) when 
they reach 7 of equation 

Re{x)^r,. (5) 

Here i\ = (1 + e)rH, while e depends on the detector 
and the 4— momentum of the particle. The presence of 
the BH event horizon is taken into by defining suitable 
boundary conditions for the kinetic distribution function 
on the hypersurface 7. For this purpose we distinguish 
between incoming and outgoing distributions on 7, /5q(x) 
and /5q(x) corresponding respectively to tIqu" < and 
Uau" > 0, where Ua is a locally radial outward 4— vector. 
Therefore the boundary conditions on 7 are specified as 
follows 

p+(x)=p(x)<S(s-s(r))^(V^-l) (6) 
Pg(x) ^ (7) 

It follows that it is possible to represent the kinetic dis- 
tribution function in the whole space time manifold in 
the form 

PG (x) = Pg (x) + (x) (8) 

where 

p±(x) = p{x)6{s-s(t))S{^^-1)x 

xe±(i?,(x) -re(s(T))) (9) 
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with 8^ respectively denoting the strong and the weak 
Heaviside functions 



and 



e+(a) 



e-(a) = 



1 for a > 
for a <0. 



1 for a > 
for a< 0. 



(10) 



(11) 



We stress that in the above boundary conditions no de- 
tailed physical model is actually introduced for the par- 
ticle loss mechanism, since all particles are assumed to 
be captured on the same hypersurface 7, independent of 
their mass, charge and state. This provides a classical 
loss model for the BH event horizon. 

Let us now consider the evolution of the the kinetic dis- 
tribution function pg(x) in external domains, i.e. outside 
the event horizon. Assuming that binary collisions are 
negligible, or can be described by means of a mean field, 
and provided that the phase space volume element is con- 
served, it follows the collisionless Boltzmann equation, or 
the Vlasov equation in the case of charged particles [1^ , 



ds 
d7 



dr'^ dp{x) dpfj dp{x) 



ds drf^ ds dp^ 



= 



(12) 



with summation understood over repeated indexes, while 



p(x) denotes pci'^) evaluated at 



r"(s(r)) and 



Po 



dro{s) 
ds 



s{t) 



This equation resumes the conser- 



vation of the probability in the relativistic phase space in 
the domain external to the event horizon. Invoking the 
Hamiltonian dynamics for the system of point particles, 
the kinetic equation takes the conservative form 



ds 
rf7 



[p(x),i/], = 0. 



(13) 



Let us now introduce the appropriate definition of kinetic 
entropy S{p) in the context of relativistic kinetic theory. 
We intend to prove that in the presence of the BH event 
horizon it satisfies an H theorem. The concept of entropy 
in relativistic kinetic theory can be formulated by direct 
extension of customary definition given in nonrelativistic 
setting [1^ fl^I [13. Il3 | . For this purpose we introduce 
the notion of kinetic entropy measured with respect to 
an observer endowed with proper time t as follows 



S{p) = ~PJ dx{s)Sis ^ s(T))<5(y^- l)p(x)lnp(x) 

^ (14) 

where p(x) is strictly positive and, in the 8-dimensional 
integral, the state vector x is parameterized with respect 
to s, with s denoting an arbitrary arc length. Here P is 
the principal value of the integral introduce in order to 
exclude from the integration domain the subset in which 



the distribution function vanishes. Hence S{p) can also 
be written: 

S{p) = -pJ dx{s)5{s - s(r))pi(x) Inp(x), (15) 
r 

where pi (x) now reads 



pi(x) - e(r(s) - r,{s))S{^^ -l)p{x{s)). (16) 

Differentiating with respect to r and introducing the in- 
variant volume element d'^rd'^p, the entropy production 
rate results manifestly proportional to the area A of the 
event horizon 



dSjp) 
dr 



-pJd^rd^pFr.rAHr-re)p\np]. (17) 



Indeed, the r.h.s represents the entropy flux across the 
event horizon. Moreover here, F" is the subdomain of 
phase space corresponding to the particle falling into the 
BH and Fj-r, is the characteristic integrating factor 



p , 



_ ds{T) / dr dr, 



dr V ds ds 



(18) 



We can write the above expression in terms of the ki- 
netic probability density evaluated at the hypersurface 
■y/u^^u^^— T, defined as /(x) = p/N. It follows 

p\np = NflnNf. (19) 

At this point we adopt a customary procedure in statis- 
tical mechanics (16.j invoking the inequality 



Nf In TV/ > TV/ - 1 



(20) 



and notice that in the subdomain of phase space F^ in 
which Frr^ > there results by definition p = 0. Hence 
it follows that in F^ 



Frr. < 



(21) 



where by construction ^^^-^ > 0. This result holds in- 



dependent of the value of ^ . Next let us introduce the 
bounded subset C F^ such that Nf results infinitesi- 
mal (of order 6 ) in the complementary set F~ \ f2 and 
moreover the ordering estimate 



P J d^rd^p\FrrJS{r-r^) 

r-\n 



NflnNf 



0{S\ii6) 



is required to hold. Manifestly this domain includes the 
set of improper points of F". In the set fi, Nf is by 
assumption positive and such that N f > 6. Therefore 
there results 

0<PJ d^rd^P \Frr, \6{r-r,)= Ms (22) 
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where Ms is a suitable finite constant. Thus one obtains 

^^>pfdhd^p\Frr \Sir-r,)[Nf -l] + 0{S\ii5) 
dr J 
n 

(23) 

The first term of the r.h.s of (|23|) can be interpreted in 
terms of the effective radial velocity of incoming particles 



Bekenstein entropy H17|l . Finally interesting features of 
the derivation are that the entropy production rate re- 
sults proportional to the area of the event horizon and 
that the formalism is independent of the detailed model 
adopted for the BH. In particular also the possible pres- 
ence of an imploding star (contracting event horizon) is 
permitted. 



V,'" = — [d^p\FrrJS{r-r,)f, (24) 
no J 
n 

while riQ is the surface number density of the incoming 
particle distribution function 

J d^rd^pfiK)Sir - ro(s(t)) = (25) 
h 

Finally we invoke the majorization 

^ > ^^iVinfjl d'rnoV^'ff^ 

-Ms + OiS\nS) (26) 

and impose that iV ^ 1 be sufficiently large to satisfy 
the inequality S > 0. We stress that inf {/ d^rnoV^^^} 
can be assumed strictly positive for non isolated BHs 
surrounded by matter. This proves the relativistic H- 
theorem 

^ > 0. (27) 
dr 

Let us briefly analyze the basic implications of our re- 
sults. First we notice that the H-theorem here obtained 
appears of general validity even if achieved in the classical 
framework and under the customary requirement iV ^ 1 
(large classical system) and in validity of the subsidiary 
condition S > 0. Indeed the result applies to BH hav- 
ing, in principle, arbitrary shape of the event horizon. 
The description adopted is purely classical both for the 
falling particles (charged or neutral 0,0,0,1131) and 
for the gravitational field and is based on the relativis- 
tic collisionless Boltzmann equation and/or the Vlasov 
equation respectively for neutral and charged particles. 
A key aspect of our formalism is the definition of suitable 
boundary conditions for the kinetic distribution function 
in order to take into account the presence of the event 
horizon. Second, the expressions for the entropy and 
entropy production rate, respectively Eas. lfni) and ((T7|l . 
can be used to determine the Bekenstein entropy for clas- 
sical BH ^ and the related generalized second law ||2J), 
although we stress that the present results are indepen- 
dent of the assumptions involved in the definition of the 
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